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We present an isotropic spin wave (ISW) theory of short-range order in Heisenberg magnets, and 
apply it to square lattice S = 1/2 and 5 = 1 antiferromagnets. Our theory has three identical 
(isotropic) spin wave modes, whereas the conventional spin wave theory has two transverse and one 
longitudinal mode. We calculate temperature dependences of various thermodynamic observables 
analytically and find good (several per cent) agreement with independently obtained numerical 
results in a broad temperature range. 



Introduction. We study short range magnetic order 
in Heisenberg antiferromagnets described by the follow- 
ing Hamiltonian: 



H — J ^ ^ SfSr 



(1) 



{rr'> 



where (rr') denotes nearest neighbor pairs. The criti- 
cal behavior of square lattice Heisenberg antiferromag- 
nets is known: they have staggered magnetic order at 
T = 0, and exponentially divergent correlation length in 
the limit T ^ 0. The purpose of our paper is not to add 
anything to understanding their asymptotic properties 
for T — > 0, but rather to develop an accurate analyt- 
ical theory for calculating observables at nonzero tem- 
perature, where the correlation length is finite and the 
system has short range, but not long range, order. Re- 
cent advances in numerical methods for quantum mag- 
nets provided us with detailed numerical data to test our 
analytical predictions. 

Our theory is applicable in one and two dimensions. 
In this paper, we focus on S =1/2 and 5*= 1 square lat- 
tice antiferromagnets, such as La2Cu04 and Sr2Cu02Cl2 
(S=l/2) and La2Ni04 (S=l). 

We obtained numerical results for using high temper- 
ature series expansions for the following quantities: (i) 
static susceptibility x(q) — M(q)/H(q) for arbitrary q, 
(ii) equal time spin correlator S'(q) = (SqS_q) for arbi- 
trary q, (iii) correlation length ^(T), (iv) uniform sus- 
ceptibility Xo{T), and (v) internal energy U{T). The lat- 
ter three quantities are derived from the former two, and 
agree with the earlier Monte Carlo [§-D and series expan- 
sions data, and with the experiment |6|H]. Since it is 
unlikely for different experimental and numerical meth- 
ods to have identical systematic errors, we believe that 
the numerical results are accurate and provide a reliable 
test of the analytical theory. 



Isotropic versus Conventional Spin- Wave The- 
ory. Magnetic excitations of the Heisenberg model are 
interacting spin waves. Their interaction (mode cou- 
pling) is essential not only for dissipation, but also for 
dynamically generating one longitudinal and two trans- 
verse modes in the low-temperature limit, as required 
by the symmetry of the ordered phase. Mode coupling 
is accurately captured by the quantum nonlinear sigma 
(QNLcr) model, but only at the expense of limiting its 
applicability to the temperature range where ^ ^ 1 (f is 
the correlation length). Here and in what follows we as- 
sume the units where lattice spacing a = 1 and exchange 
constant J = 1. 

Recent detailed calculations ^ show that for T > ps 
{ps ^ JS^ is spin stiffness), where most of the experimen- 
tal and numerical data exists, the effect of mode coupling 
is rather weak. This observation led us to explore approx- 
imations which ignore mode coupling, in which case the 
theory no longer requires that ^ 3> 1. 

Spin Wave Spectrum. We linearize the equations of 
motion written in terms of variables defined on nearest 
neighbor bonds: 



Sr' , 



rtlrr' = S 
rirr' = Sr — Sr' , 
Ir 



(2) 



which form an orthogonal triad for each bond, and seek 
canonical variables as linear combinations of (urr' , Irr' ) • 
According to the exact solution of the classical Heisen- 
berg chain by Fisher JlO| , the thermodynamic properties 
of a chain are identical to those of an isolated pair of 
classical spins. Because (nrr',lrr') are exact canonical 
variables for an isolated classical spin pair, a linear the- 
ory in terms of these variables must be able to reproduce 
exactly the thermodynamics of a classical chain at arbi- 
trary temperature. Furthermore, for a system which is 
sufficiently close to ferro- or antiferromagnetic long range 



1 



order (i.e. ^ 1), the variable set (||) is suitable for the 
derivation of the corresponding ferro- or antiferromag- 
netic quantum nonlinear sigma model, which makes the 
correspondence between the sigma model theory and our 
isotropic spin wave theory more transparent. 

The equations of motion written for iirr' and Irr' on 
nearest neighbor bonds 



dt 

dlrr' 

dt 



= ^ ^ (lr,r+5 — Ir'.r'+iS) 

S 



(3) 



where Crr' = SrSr', are linearized by replacing 



Cr+S,r' 

s 



ZCl, Cr 



C2, 



(4) 



in order to have (rirr' , Irr' ) as canonical variables of the 
resulting linear theory. Here Z is the coordination num- 
ber, and sums are taken over all nearest neighbors. Ci 
and C2 do not depend on the bond rr' because they must 
obey the translational symmetry of the paramagnetic 
phase, and are determined self-consistently, in a manner 
similar to calculating the mass term in 1 /N expansion of 
the quantum nonlinear sigma model; they are not neces- 
sarily equal to the equal time averages of the spin oper- 
ators in Eqs.(|4|). The remaining fluctuation component 
becomes the mode coupling term written as |np ~ 1 
in the quantum nonlinear sigma model. In our linear 
spin wave theory, the fluctuations of Crr' around self- 
consistent averages Ci and C2 are ignored, which leads 
to the following spectrum of nonintcracting spin waves: 



Z. 



y(l-7,)(l + 07g), 



where Z^ ^ Z^Ci, 9 = -C2/C1, and 



7g = ^Xl^^P(^'l^)- 



(5) 



(6) 



For the square lattice, 7g = ^{cos{qx) + coa{qy)). 

Our method of deriving the spin wave spectrum is sim- 
ilar in spirit to the work of Villain and Haldane [|2| , 
but uses an expressly isotropic set of variables. Starykh 
recently pointed out an elegant alternative method of 
derivation, based on frequency moments of the dynami- 
cal susceptibility, which will be described elsewhere p3| . 
Villain Young and Shastry and others earlier 
obtained the same as a function of wavevector, but 
their Z^{T) and 0{T) are different from ours. 

A similar form of e„ arises in the equations of motion 



closure methods (also called decoupling methods |15|) 
that replace A^r' , -Brr' in Eq- (|^) by their equal time av- 
erages. These methods are not equivalent to our theory, 



e.g. they do not reproduce temperature dependences pre- 
dicted by the exact solution for 5* = 00. 

The dynamical spin susceptibility is derived using the 
standard quantization procedure: 



/ , 4 1 - 7t7 

x(q,«^«) = -i^ux -2—2 



n ' q 



where 



(7) 



(8) 



is the internal energy per spin. x(q, iw„) is subject to 
the following two constraints: 

/ J^xic^,^u;n)^\siS + l), (9) 



(277)° 



(2^) 



(10) 



where ujn = 2'KnT is the set of Matsubara frequencies, 
and Eq.(|l^) follows from Eq.(|). 

In the classical (S* = 00) limit, Eqs.(p[jlO|) reproduce 
the exact solution of the classical Heisenberg chain by 
Fisher [0; they remain exact for Bethe lattices with 
arbitrary coordination number. 

In two dimensions, Eqs.(|^|o|) depend on q only 
through 7q, a prediction which is not exact for lattices 
which contain loops, such as the square lattice. We have 
verifled that numerically calculated observables for the 
Heisenberg models, for example x(q) and S'(q), depend 
on q only through 7^ to several per cent accuracy, which 
is consistent with the small probability, 1/32 « 3%, for 
a four-step path on the square lattice to form a loop. 

Analytical Calculation of the Internal Energy. 
The closure of Eqs.(|-0) requires one additional con- 
straint on the three variables Z^{T), 9{T), and U{T). In 
his exact solution. Fisher uses a special property of the 
S = 00 model, namely that the internal energy U (T) for 
a chain is the same as for an isolated spin pair. Because 
this method does not apply for finite spin, we use a dif- 
ferent method based on the same assumption of linearity 
(no mode coupling) as we used to derive eq. 

If a linear spin wave theory is a good approximation, 
the internal energy at temperatures where occupation 
numbers are sufficiently small must be given by (Ander- 
son [0): 



(27r)^ exp(eO/r) - l' 



UiT) = Uo+N 
where the ground state energy Uq is ]16|-|1^]: 

.„.4...(.,i_£.,<l_^,o(± 



(11) 



(12) 
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the r = spin wave spectrum: 
I -Id 



T" S(q) and x(q) vs. q 



2S 



O 



and 



Id = 



(27r)" 



1 

2Z' 



(13) 



(14) 



The asymptotic T ^ spin wave theory [|l^,^, and 
hence the quantum nonhnear sigma model theory for the 
renormahzed classical regime predict two gapless 

transverse spin wave modes per one magnetic unit cell, 
or per two spins, i.e. TV = 1 in Eq.(ll). 

Our theory has no mode coupling, and therefore 
does not distinguish between longitudinal and transverse 
modes (such an approximation is valid for quantum spin 
models, but not for classical models where iV = 1 should 
be used); therefore, it must have three equivalent har- 
monic oscillator modes per two spins, or TV = 3/2 per 
spin in Eq.(|ll|). This can be seen, e.g., from the follow- 
ing argument: if each of the spins is assigned to a dimer 
which it shares with one of its neighbors, one set of canon- 
ical variables (rirr' , Irr' ) per each dimer correctly counts 
the degrees of freedom, three per dimer, or = 3/2 per 
unit cell. 

The crossover between = 1 and N = 3/2 regimes is 
expected to occur when the temperature is of order lon- 
gitudinal mode gap. Figure ^ shows that our prediction 
of = 3/2 is consistent with the data down to the low- 
est temperatures where numerical results are available, 
T = 0.25J for S = 1/2 and T = J for S = 1. 
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FIG. 1. Numerical results for the internal energy mea- 
sured with respect to the ground state ?7(r) — Uq (squares, 
Ref. j^; circles, this work) are compared to Eqs.( pl[]l4[ ) with- 
out adjustable parameters. Solid line corresponds to A'' = 3/2 
and dashed line to A'^ = 1 expected in the renormalized clas- 
sical regime. 
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FIG. 2. Numerical series expansions results for the equal 
time correlator r~^5'(q) (circles) and the static susceptibility 
x(q) (diamonds) across the Brillouin zone are compared to 
our analytical theory Eqs.(jl^,|l^ without adjustable param- 
eters. 
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FIG. 3. Numerical data for the correlation length ^, 
bulk susceptibility xoi and normalized Lorentzian arnplitude 
Sq/iTf) for 5" = 1/2 and 5" = 1 models, from Ref. ||| (dia- 
monds), Ref. (squares), and this work (circles). Solid Hues 
is the theory presented here; dashed lines is renormalized clas- 
sical theory [|l^ with numerical constants calculated in [po[|9[ . 
There are no adjustable parameters; note however that our 
theoretical curves are completely analytical below the temper- 
ature shown by arrow, and use numerically calculated U{T) 
as input above that temperature. 
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Comparison With the Numerical Data. 

static susceptibility 



X(q,0) = -^ 



1 



3Z2 1 + 9jq 

has exactly the same q-dependence as in the exact solu- 
tion of the classical Heisenberg chain, and is Lorentzian 
near the ordering wavevector Q. Accordingly, the critical 
dimension of spin in this theory is equal to its mean-field 
value 7] — 0. These predictions compare well with the 
actual nearly Lorentzian shape of x(q, a; = 0) and a very 
small value of 77 = 0.023 known from numerical cal- 
culations. In contrast with the static susceptibility, the 
prediction for the equal-time correlator 



S{q) 



3 eg tanh(e,/2r) 



The cases, the agreement is within few per cent in the range of 
applicability, which for these two models approximately 
corresponds to the correlation length of less than ten lat- 

(15) tice spacings. We expect this theory to be applicable for 
a variety of other magnets, and are planning to pursue 
further studies in this direction. 
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does not have Lorentzian shape near Q. These analyti- 
cally calculated x(q) and S'(q) are plotted in Fig.^ as a 
function of wavevector for the lowest temperature where 
numerical results for q-dependence are available for com- 
parison. We find good agreement of analytical and nu- 
merical curves, achieved without adjustable parameters. 

Analytical and numerical results for temperature de- 
pendences of the correlation length ^(T) = ^^9/{l — 9) 
(defined from the long-distance decay of correlations 
~ exp(— i?/^)), bulk susceptibility xo{T), and normal- 
ized Lorentzian amplitude Sq / (T^^) are plotted in Fig||. 
Good agreement between numerical and analytical re- 
sults is observed in most of the temperature range, except 
for the lowest temperatures for 5* = 1/2 (roughly ^ > 10) 
where all three quantities begin to deviate from the theo- 
retical predictions. The deviations are most likely caused 
by a crossover to the renormalized classical behavior. 
Evidently, the crossover is completed at still lower tem- 
peratures, where no numerical data is currently available; 
analysis of the experiment remains the only currently 
available alternative for studying this crossover further. 

Summary. We have developed an isotropic spin 
wave theory of short range magnetic order and compared 
it with the independently obtained numerical data for 
5 = 1/2 and S = 1 square lattice antiferromagnets. Un- 
like the quantum nonlinear sigma or any other continuous 
model, our theory does not require the correlation length 
to be much longer than the lattice spacing. Furthermore, 
because our theory is valid for all wavevectors through- 
out the Brillouin zone, it is suitable for calculating those 
quantities that are not dominated by correlations near 
the ordering wavevector. Its primary drawback, com- 
pared to the quantum nonlinear sigma model, is linear- 
ity and the resulting lack of mode coupling or dissipation 
(the continuous limit of our theory is the 0(oo) quan- 
tum nonlinear sigma model, and not the physical 0(3) 
model) . 

We found that the analytical theory agrees with the 
numerical data without adjustable parameters. In most 
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